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1.  Introduction 

To  defeat  an  incoming  projectile,  an  active  tamor  protection  system  must  be  able 
to  estimate  the  trajectory  of  the  incoming  projectile.  Many  sensor  systems  have  been 
developed  using  angle  of  arrival  (AoA)  information  from  two  separate  receivers.  As  a 
rule  of  thumb  the  error  associated  with  these  systems  is  a  function  of  range  to  the  tar¬ 
get  over  the  separation  of  the  two  receivers.  As  this  ratio  becomes  large,  the  estimate 
loses  its  value.  On  a  tank  the  maximum  separation  distance  between  receivers  is  less 
than  four  meters.  An  AoA  system  installed  on  a  tank  would  begin  to  lose  its  effective¬ 
ness  at  ranges  in  excess  of  eight  meters.  Obviously,  eight  meters  does  not  allow 
enough  reaction  time  for  an  active  protection  system  to  respond. 

H.  Bruce  Wallace  of  the  Ballistic  Research  Laboratory  (BRL)  has  proposed  that 
range  information  also  be  utilized  in  order  to  have  a  worthwhile  target  location  esti¬ 
mate.  There  is  a  large  class  of  p>otential  sensors  fitting  this  description.  These  can  be 
classTied  according  to  the  type  of  electronic  processing  they  use  and  the  parameters 
that  are  important  for  that  type  of  processing.  The  goal  here  is  to  find  a  general  model 
that  allows  for  various  levels  of  detail  in  the  performance  evaluation  of  this  class  of 
sensor  systems. 

This  paper  develops  a  performance  model  for  a  system  using  AoA  and  range 
information  to  estimate  the  location  of  a  target.  The  ideas  used  to  derive  and  validate 
the  performance  model  are  presented,  and  then  a  procedure  for  evaluating  specific 
sensor  systems  is  discussed. 

First  the  two  dimensional  case  is  examined  in  detail.  Then  the  results  of  applying 
the  same  ideas  to  the  three  dimensional  case  are  presented.  Both  models  were  vali¬ 
dated  through  the  use  of  a  simulation.  The  procedure  used  to  validate  the  three 
dimensional  model  is  presented.  The  next  two  sections  present  some  electronic 
models  from  the  literature  associated  with  the  accuracy  of  angle  of  arrival  and  range 
measurements.  Finally,  a  method  for  the  evaluation  of  a  particular  system  is  given. 

2.  XY  Covariance  Model 

In  tliis  section  a  performance  model  for  a  system  processing  one  range  measure¬ 
ment  and  one  AoA  measurement  is  presented.  The  meacurementr  qre  from  a  polar 
coordinate  system  but  the  system  performance  is  to  be  in  the  Cartesian  coordinate 
system.  The  focus  of  the  discussion  is  on  how  the  measurement  errors  effect  target 
estimation  in  the  Cartesian  plane.  The  relationship  between  an  xy  location  and  an 
AoA  range  coordinate  is  straight  forward. 

x  =  R  cos  ^  +  Xp  (1) 

y  =  R  sin  0  y^ 

when  R  is  the  range  to  the  target  and  6  is  the  AoA  to  the  target.  Without  any  loss  of 
generality,  we  assume  that  x^=0  and  yo=0 . 

A  common  approach  for  relating  measurement  errors  to  the  xy  domain  is 
through  the  partial  derivatives  of  the  location  with  respect  to  the  measured  quantities 
(Reference  1  and  2).  In  the  x  direction  the  changes  caused  by  range  measurement  are 
AR  cos  0.  Changes  caused  by  increasing  the  angle  move  x  closer  to  the  origin  by  the 
amc  ant  R  sin  (Atf)  sin(^).  In  order  io  derive  the  model,  several  assumptions  will  be 
made.  Later  it  will  be  shown  that  the  assumptions  are  statistically  reasonable. 
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Perturbations  in  x  due  to  measurement  errors  are  described  by 

X  +  A-x  =  (R  +  AR)  cos  {6  +  A0). 

After  expanding  the  cosine  term  we  have: 

(R  +  AR)  cos  (6  +  A^)  =  R  cos  6  cos  Afl  -  R  sin  0  sin  + 

AR  cos  $  cos  A0  -  AR  sin  0  sin  A0. 

Since  A0  is  small  we  will  assume  cos  A0  =  1  ;  also  since  A0  is  small,  we  assume 
sin  A0  =  A0  -,  and  finally  we  assume  A0  AR  =  0 .  So  we  now  have: 

(R  +  AR)  cos  (0  +  A0)  =  R  cos  0  +  AR  cos  0  -  RA0  sin  0. 

Recalling  that  x  =  Rcos5  we  have 

Ax  =  AR  cos  0  -  R  A0  sin  0. 

By  a  similar  argument  it  can  be  shown  that 

Ay  =  A  R  sin  0  +  RA0  cos  0. 

Note  that  taking  partials  of  both  the  x  and  y  values  in  Equation  1  results  in  the  above 
set  of  equations.  TEius,  the  terms  w'e  have  ignored  correspond  to  higher  order  differen¬ 
tials.  We  will  assume  the  measurement  errors  have  the  following  properties 

E  (AR)  =  0 

E  iA0)  =  0 

E  (AR^)  =  (2) 

E  (A^)  = 

E  (AR  A0)  =  0. 

Perturbations  of  x  and  y,  can  be  expressed  in  terms  of  measured  quantities  as : 

Ax  =  AR  cos  0  -  A0  R  sin  0 
Ay  =  AR  sin  0  -t-  A0  R  cos  0 

(Ax)^  =  (AR)^  cos^  0  -t-  {A0)^  R^  sin^  0  -  2  AR  A0  R  sin  0  cos0  (3) 

(Ay)^  =  (AR)^  sin^  0  +  (A0)^  R^  cos^  ^  -f-  2  AR  A0  R  sin  0  cos  0 

Ax  Ay  =  AR^  sin  0  cos  0  -  A^  R^  sin  0  cos  0  -  AR  A0  R  sir?  0  +  AR  A0  R  cos^  0  . 

Taking  the  expectations  of  Equation  3  and  combining  this  with  Equations  2,  we  get  the 
following  expressions  for  the  variance  and  covariance  of  the  target  location. 

E  (Ax)  =  0 

E  (Ay)  =  0 

E  (Ax^)  =  cos^  0  -t-  R^  sin^  0  (4) 

E  (Ay^)  =  sin‘  0  +  R^  cos^  0 

E  (Ax  Ay)  =  .5  sin  2  0  -  o^g  R^  sin  2 


These  equations  describe  the  error  matrix  for  the  system  in  terms  of  the  xy  coordinate 
svGiem.  In  matrix  notation,  the  covariance  matrix  of  the  estimated  position  would  be 


E  (Ax^)  E  (Ax  Ay) 

A 

E  (Ax  Ay)  E  (Ay^)  - 

.'“’A  . 

By  rotating  the  coordinates  through  the  angle  a  defined  by 

2  E  (Ax  Ay) 

tan  2  a  =  - - - 

E  (Ax^)  -  E  (Ay^) 

we  decouple  the  system  and  find  the  major  and  minor  axis.  (This  is  derived  in  Appen¬ 
dix  A). 

From  Appendix  A  we  ha\  e 

Mair-r  Axis  =  —  |e  (Ax“)  E  (Ay^)  +  ((E  (Ax^)  -  E  (Ay)^)^  +  4  E  (Ax  Ay)*) 

(5) 

Minor  Axis  =  | E  (Ax^)  +  E  (Ay^)  -  ( ( E  (Ax^)  -  E  (Ay)^)^  -f-  4  E  (Ax  Ay)*) ' 


TTie  following  equalities  were  used  to  rewrite  Equation  5. 

E  (Ax*)  -I-  E  (Ay*)  =  a*j^  sin^  6  +  R*  cos^  6  -i-  cos^  6  +  R*  sin*  6 
=  cr*j^  (cos*  6  -I-  sin*  $)  +  o'g  R*  (cos^  B  +  sin^  6) 

=  o\  +  o  gR^ 

E  (Ax*)  -  E  (Ay*)  =  (cos^  6  -  sin*  6)  +  o^g  R^  (sin^  6  -  cos^  6) 

=  cos  26-o^g  R^  cos  2 0 
E(Ax  Ay)  =  .5  (  sin  2  0  -  a^g  R^  sin  2  ^ ) 

Thus  the  major  axis  and  minor  axis  are  defined  by  : 

^  ( ( a*R  R*  a*^)  ±  ( (-a*R  -H  o\  R^)^  cos^  2  6  ^  (a\  -  R^)^  sin^  2  0)  ) 
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The  final  simplification  of  the  above  equation  results  in  principal  components  of 
and  R  CT  j .  In  the  case  of  uncorrelated  measurement  errors,  the  axes  describing 
the  covariance  structure  are  directly  related  to  the  measurement  errors. 

For  the  two  dimensional  case,  performance  can  be  discussed  in  terms  of  x  and  y 
or  in  terms  of  principal  components.  Although  the  ideas  are  the  same  in  three 


dimensions  the  relationships  are  more  complex  and  it  is  difficult  to  give  a  purely  tri¬ 
gonometric  explanation. 

3.  XYZ  Covariance  Model 

In  three  dimensions  a  covariance  model  can  be  derived  from  range,  azimuth 
angle,  and  elevation  measurements.  The  elevation  angle,  0 ,  is  measured  from  the  posi¬ 
tive  z  axis  and  the  azimuth  angle,  <^,  is  measured  from  the  x  axis  toward  the  positive  y 
axis.  Note  that  conventional  notation  is  being  used  and  6  has  a  different  meaning  in 
this  section  than  its  previous  meaning.  Using  the  same  simplification  arguments  as  for 
the  XY  Covariance  Model  the  following  set  of  measurement  error  transformation 
equations  is  obtained. 

Ax  =  AR  sin  6  cos  4>  +  A^  cos  6  cos  <^  -  A(^  R  sin  0  sin  <i> 

Ay  =  AR  sin  0  sin  ip  A0R  cos  ^  sin  +  A^  R  sin  0  cos  4> 

Az  =  AR  cos  0  -  A0  Rsin  0 


From  these  equations  the  covariance  model  can  be  found,  the  quantities  of  interest 
are  as  follows. 

cr  ^  =  (7  sin  0  cos"  (^  +  cr  ^  R"  cos"  0  cos"  <;i>  -t-  a  ^  R"  sin  0  sin"  4> 

sin  0  sin"  p  -t-  R“  cos  0  sin"  -i-  a  ^  R"  sin  0  cos  <p 

2  ^  2  2 
cr  7  =  cr"j^  cos"  ^  -I-  (7  ^  R"  sin  0 

^XY  ~  ^  ^  ^  “*■  *^^9  ^  ^  ^ '  0sin4>  cos  4> 

2  O'*) 

cTxz  =  <7  p  sin  6)  cos  0  cos  (p  -  o'g  R"  sin  0  cos  0  cos  <t> 

'y  2  . 

rjyz  =  ct"j^  sin  0  cos  0  sin  (;;ii  -  a  ^  R"  sin  0  cos  0  sin 

By  finding  the  eigenvalues  and  eigenvectors  associated  with  this  matrix,  the  principal 
components  and  orientation  can  be  found. 

4.  Noncolocated  Sensors 


The  detection  elements  of  a  sensor  system  are  not  usually  located  at  the  same 
position.  When  the  target  is  far  away  the  errors  associated  with  assuming  colocation 
are  negligible.  For  an  active  protection  system  the  distance  to  the  target  will  become 
small,  hopiefully  not  too  small,  and  problems  will  result  from  assuming  colocation.  In 
this  section,  the  effects  of  separating  the  range  and  AoA  detectors  will  be  discussed. 
The  rsv'O  topics  of  this  section  are  an  estimation  technique  and  a  general  method  for 
finding  the  covariance  structure  around  a  target  location. 


Let  (x,y,z)  be  the  position  of  the  range  sensor. 


Let  (0,0,0)  be  the  piosition  of  the  interferometer. 


The  distance  from  the  interferometer  to  the  target,  K,  has  the  following  xyz  com¬ 
ponents 

(K  sin5cos<?!i,  Ksin^sin<^,  Kcosd) 
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The  sum  of  the  squares  in  the  x,  y,  and  z  directions  should  equal  the  square  of  the 
measured  range.  We  proceed  by  finding  K  in  terms  of  R. 


R^  =  (Ksin^  cos4>  -  x)^+  (Ksin^  simp  -  y)^+(Kcostf  -  z)^ 

=K^sin^^  cosV  *  2Ksii^  cos<^  x  +  x^ 

+  K^sin"0  sinV  -  2Ksintf  siiK^  y  + 

+  K'cos^^  -  2Kcos^'  z  +  z^ 

2 

In  this  case  the  coefficient  for  K  is  one.  K  is  found  by  using  the  quadratic  formula. 

-b±(b^-  4c)'^ 

K= - ^ - — 

2 

where  o=x  +y  +z  -  R" 

and  b  =  -2zcos/9  -  2xsin0  cos(p  -  2ysin0  sincP . 

The  next  goal  is  to  find  a  general  method  for  determining  the  covariance  struc¬ 
ture  around  a  target  location.  The  measurement  error  structure  can  be  determined  a 
priori  from  performance  knowledge  of  the  particular  sensors.  The  typical  case  is 
uncorrelated  measurements  between  the  sensors;  and  leads  to  a  diagonal  matrix 
describing  the  error  structure  in  the  measurement  space.  The  problem  is  to  find  a 
representation  of  this  error  structure  in  a  different  reference  system.  The  associated 
theory  is  discussed  by  Dempster  (3).  The  procedure  is  to  find  representations  of  the 
reference  basis  in  terms  of  the  measurement  error  basis.  The  matrix  P  is  formed  by 
concatenating  these  column  vectors.  A  vector  in  the  reference  frame,  R,  has  the 
representation  PR  in  the  measurement  frame.  To  find  the  covariance  structure  in  the 
reference  frame  consider  any  vectors,  Rj  and  R2,  from  this  space.  To  find  the  correla¬ 
tion  between  the  two  vectors,  the  inner  product  is  taken;  however,  the  inner  product  is 
defined  in  the  measurement  space  by  the  measurement  covariance  matrix,  E.  The  fol¬ 
lowing  equation  illustrates  the  process. 

<Rj,R2>  =  <PRj,PR2> 

=  (PRj)'EPR2 
=  R'jP‘EPR2 

From  this,  the  covariance  in  the  reference  space  can  be  seen  to  be  P’EP.  Unfor¬ 
tunately,  in  the  typical  case  it  is  not  possible  to  find  P  directly.  The  matrix  Q  (  =  P’^) 
can  be  found  by  finding  the  representations  of  the  measurement  basis  in  terms  of  the 
reference  basis,  and  concatenating  these  as  columns.  By  finding  the  inverse  of  Q,  the 
desired  matrix  P  is  found.  Notice  if  Q  is  ill  conditioned,  P  will  include  some  laige 
values.  If  P  is  an  orthonormal  set,  then  the  process  is  a  rotation  and  P*  is  unified  with 
P‘^  =  O. 

The  process  will  be  illustrated  for  the  previously  discussed  three  dimensional 
case.  The  errors  associated  with  each  measurement  will  be  described  in  sequence. 
Errors  in  the  azimuth  measurement,  <P,  result  in  perturbations  in  the  xy  plane 
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perpendicular  to  the  line  defined  by  4>  .  ITie  normalized  vector  in  this  direction  is 
(-cos</> ,  sin<^ ,  0 ).  The  vector  associated  with  the  elevation  measurement,  B,  is  perpen¬ 
dicular  to  the  ray  going  through  the  target  and  contained  in  the  plane  defined  by  the  xy 
direction  to  the  target  with  z  unrestricted.  Recall  that  B  is  measured  from  the  positive  z 
axis  toward  the  xy  plane.  The  normalized  vector  is  ( -costf  sin*^ ,  -costf  cos<^ ,  sinB ).  The 
direction  of  the  ran<^e  measurement  error  is  along  the  ray  cormecting  the  range  sensor 
with  the  target  and  is  (  sin^  sin*^ ,  sin^  cos^ ,  cos^ ).  It  is  straightforward  to  verify  that 
these  form  an  orthonormal  set.  Concatenating  these  vectors  as  columns  forms  the 
matrix  Q,  or  in  this  case  P*.  The  magnitudes  of  the  measurement  errors  need  to  be 
placed  on  the  diagonal  of  a  matrix.  The  variance  of  the  error  associated  with  (f)  is 
<T^^(Rsin5)^,  where  R  is  the  distance  to  the  target.  The  range  error’s  magnitude  is 
Finally,  the  variance  associated  with  errors  in  measuring  B  is  (Ra^^  .  It  can  be  verified 
that  this  gives  the  same  result  as  Section  3  when  the  matrices  P  E  P  are  multiplied. 
This  approach  is  general  in  that  it  can  be  used  for  any  nondegenerate  situation. 

The  case  where  the  range  sensor  is  not  colocated  with  the  interferometers  will  be 
discussed  next.  The  interferometer  measurements  are  the  same  as  discussed  above. 
TTie  range  measurement  direction  is  found  by  normalizing  the  ray  between  the  range 
sensor  and  the  target.  If  D  is  the  distance  to  the  target  and  (xj^ ,  yj^ ,  Zj^)  are  the 
respective  x,  y,  and  z  distances,  then  the  normalized  vector  is  found  by  dividing  each  of 
these  quantities  by  D.  Note  in  this  case  the  set  of  normalized  vectors  will  rarely  form 
an  orthonormal  set  and  thus  Q’^  must  be  calculated  to  find  P.  It  is  worth  noting  that  as 
the  target  gets  closer  to  the  sensors,  the  projection  of  the  range  measurement  into  the 
space  defined  by  the  two  angle  measurements  will  increase.  This  increase  in  collinear- 
ity  becomes  pronounced  as  the  distance  to  the  target  approaches  the  separation  dis¬ 
tance  of  the  sensors. 


5.  Model  Validation  Effort 

Both  the  xy  and  the  xyz  models  were  compared  with  simulation  data  to  verify 
their  performance.  Using  a  Gaussian  random  number  generator,  errors  for  range  and 
azimuth  or  for  range,  azimuth,  and  elevation  were  generated.  These  errors  were 
added  to  the  true  values  and  then  the  estimated  position  was  calculated  from  the  cor¬ 
rupted  values.  Using  ten  thousand  such  points,  the  covariance  of  the  target  position,  S, 
was  calculated  and  then  compared  to  the  covariance  predicted  by  the  model,  E.  The 
code  designed  to  jxrform  the  simulation  is  included  as  Appendix  B. 

As  a  first  test,  the  determinants  were  compared  to  see  if  they  were  in  agreement. 
A  test  based  on  the  asymptotic  distribution  of  the  sample  covariance  was  used  for  this. 
The  assumption  under  the  null  hypothesis  is  that 

det(S) 

—^-1) 
det(E) 

is  distributed  as  N(0,2p) 


where 

p  is  the  dimension  of  the  matrix 
n  is  the  degrees  of  freedom. 


0 


In  this  case,  a  nonrejection  of  the  null  hypothesis  indicates  that  global  error  associated 
with  each  covariance  structure  is  similar;  however,  an  acceptance  of  the  model  based 
upon  a  determinant  test  is  not  conclusive.  These  tests  have  no  sensitivity  to  the  orien¬ 
tation  of  the  covariance  structure,  and  will  not  detect  certain  differences  in  the  magni¬ 
tude  of  the  major  components  of  the  matrix.  For  example,  a  diagonal  matrix  with 
components  (1  ,1)  has  the  same  determinant  as  the  diagonal  matrix  with  components 
(.1, 10)  or  the  matrix  with  row  one  elements  (5,  3)  and  row  two  elements  (3, 2). 

The  statistical  properties  of  a  covariance  structure  are  described  by  the  Wishart 
distribution.  Dempster(3)  discusses  the  related  theory.  The  use  of  this  distribution  in 
the  three  dimensional  case  leads  to  six  independent  tests,  each  one  checking  a 
separate  component  of  the  covariance  matrix.  Note  that  the  off-diagonal  terms  are 
based  on  estimates  of  the  diagonal  terms  and  thus  have  more  uncertainty  associated 
wth  them.  The  procedure  used  was  as  follows: 

1.  Calculate  the  error  structure  from  the  model. 

2.  Calculate  the  error  structure  from  the  simulated  data. 

3.  Find  the  normalizing  transformation  based  on  the  model. 

4.  Apply  the  normalizing  transformation  to  the  result  of  step  2. 

5.  Test  the  resulting  matrix  to  see  if  it  is  statistically  equivalent  to  the 

identity  matrix  (see  Dempster(3)). 

In  each  of  the  cases  investigated  the  model  and  the  simulated  data  produced  covari¬ 
ance  structures  that  were  statistically  the  same. 

6.  AoA  Errors 

The  Pvo  models  included  in  this  section  are  based  entirely  on  thermal  noise  and 
should  be  used  as  the  best  case  situations  for  AoA  errors.  When  additional  sources  of 
error  are  modeled  it  is  usually  correct  to  take  their  root  mean  square  with  the  error 
due  to  thermal  noise.  Dr.  Alexander  in  (2)  gives  the  following  two  equations  for  relat¬ 
ing  electronic  parameters  to  the  angle  measurement  error. 

For  pulsed  AoA  processing  the  thermal  noise,  ,  of  a  phase  interferometer  is 
given  by 

(360/27r)C 

o  u  =  - . 

27rfdcos^4(S/N)jNTG 


For  an  amplitude  monopulse  the  thermal  noise  is 


25.4  0^ 


‘^th  = 


1/2 


where, 

0  is  the  angle  form  boresight  to  the  target, 
c  is  the  propagation  velocity  (M/S), 
f  is  the  RF  carrier  (Hz) . 
d  is  the  spacing  between  the  receiving  antennas. 
0Q  is  the  antenna  half  power  beam  width  (deg). 
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^/^’iNTG  is  the  integrator  output  of  signal  to  noise  ratio. 


7.  Range  Errors 

Errors  in  range  depend  on  the  ability  to  measure  the  time  of  arrival  of  a  given 
pulse.  The  range  error  is  =  C/2  where  is  the  time  error.  The  time  error  is 
dependent  on  pulse  type,  the  following  are  suggested  by  Skolnik(4): 

For  a  rectangular  pulse 


4/3E/No 


where  r  is  the  pulse  width 
13  is  the  bandwidth 
E  is  the  received  signal  energy 
Nq  is  the  noise  per  unit  bandwidth. 


For  a  trapezoidal  pulse 


+  3  Tj  T, 
6E/N„ 


Tj  is  the  time  duration  of  the  pulse 

T2  is  the  rise  and  fall  time  of  the  top  of  the  pulse. 


For  a  triangular  pulse 


\  = 

/12  (2 


1.384r 

For  a  Gaussian  pulse  of  the  form  s  (t)  =  exp  (- - ) 

2 

T 


7r/?(2E/No)' 


If  the  pulse  has  the  form 


sin  (tt  ^  r) 


For  Continuous  Wave,  the  error  is 


4  7rAf(2  E/Nq)^/^ 
where  A  f  is  the  difference  between  the  two  frequencies. 
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The  specific  electronic  model  will  vary  depending  on  the  method  itsed  to  extract 
information  from  the  signal.  The  models  included  above  only  represent  some  of  the 
common  techniques.  The  important  thing  is  that  they  relate  specific  electronic  or 
geometric  parameters  to  the  measurement  errors  associated  with  either  range  or 
angle  of  arrival. 

8.  Analysis  Method 

In  using  hese  models  to  analyze  the  performance  of  a  system,  the  following  steps 
must  be  taken. 

1.  Choose  a  AoA  error  model  and  assign  values  to  the  parameters. 

2.  Choose  a  range  error  model  and  assign  values  to  the  parameters. 

3.  Use  the  range  and  AoA  errors  as  input  to  the  xyz  location  model. 

By  following  this  procedure  at  a  number  of  different  points  a  system’s  perfor¬ 
mance  can  be  presented  as  a  function  of  target  location. 

9.  Conclusion 

This  application  is  t>'pical  of  the  error  analysis  approach  used  in  many  engineer¬ 
ing  studies.  In  this  case  ,  the  extra  step  of  checking  the  statistical  validity  of  the  model 
was  included.  The  dominant  feature  of  this  approach  to  system  analysis  is  to  start  with 
the  measurement  errors  and  follow  them  as  they  propagate  through  the  system  and 
degrade  the  ideal  system  performance. 

The  models  presented  herein  can  be  used  to  evaluate  the  performance  of  many 
range  -  angle  sensor  systems.  The  performance  can  be  based  on  specific  electronic 
parameters  such  as  frequency,  pulse  shape,  or  on  more  general  specifications  such  as 
three  degree  angle  with  five  percent  range  errors.  This  work  could  be  continued  by 
designing  a  software  package  that  includes  the  selection  of  the  possible  options. 

This  model  may  be  used  to  simulate  a  stream  of  observations  or  to  provide  infor¬ 
mation  about  the  error  associated  with  an  observation.  The  performance  of  filtering 
techniques  is  enhanced  by  knowledge  of  error  associated  with  given  observations.  In 
evaluating  candidate  systems  for  an  active  armor  protection  system  this  model,  covers 
systems  processing  range  and  AoA  informatioa 


9 


Intentionally  left  blank. 
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Appendix  A 

Principle  Components  for  Two  Dimensions 
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appendix  a 

We  wish  to  find  the  angle  q  of  rotation  in  the  coordinate  system  to  decouple  the 
system. 

Ax '  =  Ax  cos  Q  +  Ay  sin  a 
Ay '  =  -Ax  sin  a  +  Ay  cos  a 


In  terms  of  the  new  coordinate  system  we  have 
E(Ax')  =  0 


E(Ay')  =  0 

E(Ax  '^)  =  E  (Ax^)  cos^  Q  +  2  E  (Ax)  (Ay)  sin  a  cos  a  +  E  (Ay^)  sin^  a 
=  E  (Ax^)  +  E  (Ay^)  -  E  (Ax^)  sin^  a  -  E  (Ay^)  cos^  a 
+  2  E  (Ax  Ay)  sin  a  cos  a  . 


Recall  that 


1 


sin  a  cos  a  =  —  sin  2  a 
2 

2  1  -  COS  2  Q 

sin^  a  - 


2 

2  1  +  cos  2  a 

cos  a  = 

2 


Then 


E  (Ax'^)  =  -  (E  (Ax^)  +  E  (Ay^))  + 
2 


1^2 
—  (E  (Ax'^)  -  E  (Ay  ))  cos  2  a 

2 


+  E  (Ax  Ay)  sin  2  a  . 


Similarly  it  can  be  shown  that 

E  (Ay'^)  =  —  (E  (Ax^)  +  E  (Ay^))  -  — (E  (Ax^)  -  E  (Ay^))  cos  2a 
2  2 

-  E  (Ax  Ay)  sin  2  a 


and 

^  2  2 
E  (Ax'  Ay')  =  -  E  (Ax  )  sin  a  cos  a  +  E  (Ax  Ay)  (cos  a  -  sin  a) 

2 

+  E  (Ay  )  sm  a  cos  a 

=  —  [E(Ay^)-E(Ax^)]  sm2a  +  E(AxAy)  cos2a. 
2 

The  covariance  term  will  be  zero  if 
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tan  2  Q  = 


2  E  (Ax  Ay) 

E  (Ax^) .  E  (Ay^) 

Note  that  for  this  angle 

2  E  (Ax  Ay) 

sin  2  a  =  - 

(4  E  (Ax  Ay)^  +  (E  (Ax^)  -  E  (Ay^))^)^/^ 

E  (Ax^)  -  E  (Ay^) 

COS  2  Q  =  - . 

(4  E  (Ax  Ayf  +  (E  (Ax^)  -  E(Ay^))^)^/^ 

Using  these  relations  it  can  be  shown  that  the  major  and  minor  axis  are  defined  by 
1/2  (E  (A\^)  +  E  (Ay^)  ±  ((E  (Ax^)  -  E(Ay"))^  +  4  E  (Ax  . 
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Appendix  B 
Validation  Program 
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APPENDIX  B 


#include  <math.h> 

#mclude  "ranvar.h" 

#mclude  "stat.h" 

main() 

{ 

float  range, r,r_sd,r_var,r_sq; 
float  theta,th,theta_sd,theta_var; 
float  psi,p,psi_sd,psi_var; 

float  x,y,z,x_sum,y_sum,z_sum,x_var,y_var,z_var; 
float  xsq,ysq.zsq,xy,xz,yz; 
float  r_sim,th_sim,psi_sim; 

float  sinth,costh,smpsi,cospsi,sin2th,cos2psi,cos2th,sin2psi; 
float  cov_xy,cov_xz,cov_yz: 

float  x2_mod,y2_mod,z2_mod,xy_mod,xz_mod,yz_niod; 
float  data_det,model_det, minus, plus, z_test,prob; 

int  i,n,seed; 

/•  initialize  the  variables  for  this  run  */ 
range  =  50; 

r_sd  =  .05*range;  /*  five  percent  range  */ 
r_var  =  r_sd*r_sd; 
r_sq=range*range; 
theta  =  M_PI/4; 

theta_sd=5*2*M_PI/360;  /’five  degree  error  elevation  */ 
theta_var  =  theta_sd*theta_sd; 
psi  =  M_PI/4; 

psi_sd=5*2*M_PI/360;  /*  five  degree  error  azimuth  */ 

psi_var =psi_sd*psi_sd; 

n  =  10000;  /*  number  of  replications  */ 

seed =23719;  /*  random  number  seed  */ 

xsq=0; 

ysq=0; 

zsq=0; 

x_sum=0; 

y_sum=0; 

z_sum=0; 

xy=0; 

xz=0; 

yz=0; 

for  (i=0;i<n;i+  +) 
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{ 

/•  rv_gaiiss  is  a  gaussian  random  number  generator 
the  first  two  sections  find  the  measured  value 
and  calculate  the  position 

cnf  is  the  cumulative  normal  distribution  function 

and  is  accurate  to  six  decimal  places 

» / 

/ 

th  =  theta +theta_sd*rvjgauss(seed); 
r = range + r_sd*rv_gauss(seed); 
p = psi + psi_sd  *  rv_gauss(seed) ; 

z  =  r*cos(th); 

x=r*sin(th)*cos(p); 

y=r*sin(th)*sin(p); 

x_sum+  =x; 
y_sum+  =y; 
z_sum+  =z; 

xsq+  =x*x; 
ysq+  =y*y; 
zsq+  =z*z; 
xy+  =x*y; 

XZ+  =x*z; 
yz+  =y*z; 

}  /•  end  of  replication  loop  */ 

x=x_sum/n; 

y=y_sum/n; 

z=z_sum/n; 

/•  the  following  values  are  the  covariance  elements  based  on 
the  simulations  data 

V 

x_var  =  (xsq-x  •x_sum) /( n- 1 ) ; 
y_var  =  (ysq-y  •y_sum)/ (n- 1); 
z_var  =  (zsq-z*z_sum)/(n- 1); 
cov_xy = (xy-x  •y_sum)/(n- 1 ) ; 
cov_xz  =  (xz-x  •  z_sum)  /  ( n- 1 ) ; 
cov_yz = (yz-y  •z_sum) /(n- 1), 

/•  the  next  section  uses  the  model  to  find  the  predicted  covariance 
structure 

V 

sinth=sin(theta); 

sin2th=sinth*sinth; 

costh=cos(theta); 
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cos2th = costh  •  costh ; 
sinpsi  =  sin(psi); 
sin2psi =sinpsi*sinpsi: 
cospsi=cos(psi); 
cos2psi =cospsi*cospsi: 


x2_mod  =  r_var*sin2th*cos2psi; 

x2_mod  +  =theta_var*r_sq*cos2th*cos2psi; 

x2_mod  +  =psi_var*r_sq*sin2th*sin2psi; 


y2_mod  =  r_var*sin2th*sin2psi  +  theta_var*r_sq*cos2th*sin2psi; 
y2_mod+  =psi_var*r_sq*sin2th*cos2psi; 

z2_niod  =  r_var  *00521  h  +  theta_var  *r_sq*siii2th: 

xy_mod  =  r_var  *sin2th*cospsi*sinpsi  +  theta_var  *r_sq*cos2th*cospsi*sinpsi; 
x}'_mod-  =psi_var*r_sq*siri2th*cospsi*sinpsi; 

xz_mod  =  r_var*sinth*costh*co5psi'theta_var*r_sq*sinth"‘costh*cospsi; 
yz_mod  =  r_var*sinth*costh*sinpsi-theta_var*r_sq*simh*costh*sinp5i; 

/•perform  determinate  test 

V 

plus  =  x_var*y_var*z_var; 
plus  +  =  cov_xy*cov_yz*cov_xz*2; 
minus  =  cov_xz*cov_xz*y_var; 
minus  +  =  cov_xy*cov_xy*z_var; 
minus  +  =  cov_yz*cov_yz*x_var; 
data_det  =  plus-minus; 

plus  =  x2_mod*y2_mod*z2_mod; 
plus  -t-  =  2*xy_mod*xz_mod*yz_mod; 
minus  =  xy_mod*xy_mod*z2_mod; 
minus  +  =  y2_mod*xz_mod*xz_mod; 
minus  +  =  x2_mod*yz_mod*yz_mod; 
model_det = plus-minus; 

z_test  =  sqrt(n-l)  *  (data  det/model  det  - 1); 
z_test  /=  sqrt(2*3); 
prob  =  cnf(z_test); 

printffO); 

printfC . . . 

printf("  Range :  %f  theta :  %f  psi :  %f', range,  theta,psi); 
printf("Od(range)  =  %fsd(theta)=  %fsd(psi)=  %f',r_sd,theta_sd,psi_sd); 
printf("Oimulated  Var(X)=  %f  model  value  was  %f'pc_varpc2_mod); 
printf("Oimulated  Var(Y)=  %f  model  value  was  %f',y_var,y2_mod); 
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printf("Oimulated  Var(Z)  =  %f  model  value  was  %f',z_var,z2_mod); 
printffOimulated  Cov(XY)  =  %f  model  value  was  %f ’,cov_xypcy_mod); 
printf("Oimulated  Cov(XZ)=  %i  model  value  was  %f',cov_xz^_mod); 
printf("Oimulated  Cov(YZ)  =  %f  model  value  was  %f',covjirz,yz_mod); 
printf("0  value  of  %f  with  probability  of  %f',z_test,prob); 

} 
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